A Bound on Binding Energies and Mass 
Renormalization in Models of Quantum 
Electrodynamics 

Elliott H. Lieb^ and Michael Loss^ 

1. Departments of Mathematics and Physics, Jadwin HaU, 
Princeton University, P. O. Box 708, Princeton, NJ 08544 

2. School of Mathematics, Georgia Tech, Atlanta, GA 30332 
March 31, 2002, 

Abstract 

We study three models of matter coupled to the ultraviolet cutoff, quantized radiation field 
and to the Coulomb potential of arbitrarily many nuclei. Two are nonrelativistic: the first uses 
the kinetic energy {p+eA{x))'^ and the second uses the Pauli-Fierz energy {p+eA{x))'^+ea-B{x). 
The third, no-pair model, is relativistic and replaces the kinetic energy with the Dirac operator 
D{A), but restricted to its positive spectral subspace, which is the "electron subspace". In each 
case we are able to give an upper bound to the binding energy - as distinct from the less difficult 
ground state energy. This implies, for the first time we believe, an estimate, albeit a crude one, 
of the mass renormalization in these theories. 

Dedicated to two masters of mathematical physics, David Ruelle and Yakov Sinai, 
on the occasion of their 65th birthdays 



1 Introduction 

There has been a great deal of recent work dedicated to the construction of a theory of ordinary bulk 

matter interacting with the quantized radiation field. In such theories the number of electrons, N, 

is usually held fixed (i.e., pair production is not allowed) and these particles interact with each 

other and with K fixed nuclei via the ordinary electrostatic Coulomb potential - in the Coulomb 

gauge. (The nuclei are fixed because they are, relative to the electrons, infinitely massive.) The 

electrons also interact with the magnetic vector potential A which is quantized and which has the 

well known quantized field energy. It is essential, however, to have an ultraviolet cutoff A in the 
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interaction of electrons and A field, for otherwise almost everything in the theory becomes infinite. 
This is not an enormous drawback since we are interested only in the low energy physics of atoms, 
molecules and photons. 

Attention has been paid mostly to the stability of matter, namely that the ground state energy 
is bounded below by C{N + K), where C is some universal constant (depending on the parameters 
of the theory but not on N and K) [||, ^, ^, 10, 13|. Almost no attention has been paid to the 
estimation of the atomic binding energy, i.e., to the difference between the ground state energies 
with and without the Coulomb potential. 

In this paper we shall consider two nonrelativistic theories and one relativistic theory. In the 
absence of quantized fields, the ground state energy (or 'self energy') of a free electron (i.e., without 
other electrons and nuclei) is zero in the nonrelativistic case and equals mc^ in the relativistic case. 
(Here, m is the unrenormalized, or "bare" mass of the elctron.) When N electrons are present, 
but without nuclei, the energy is still zero (or Nmc^) because the electrons can move infinitely far 
apart. Therefore, in the presence of nuclei, the ground state energy is equal to the binding energy 
(or equals the binding energy plus Nmc^) when there is no quantization of the A field. 

The situation changes dramatically when the A field is quantized. The self-energy of a free 
electron (i.e., the ground state energy without Coulomb potentials but with the quantized A field) 
is large if A is not too small. If the fine structure constant a = /he is not too large (e.g., 1/137) 
and if the nuclear charges are not too large then the change in the ground state energy is not very 
large. Thus, the binding energy is the difference of two large quantities and its calculation is like 
"looking for a needle in a haystack" . 

All three models use 'minimal coupling', i.e., p is replaced by eA{x)/c in the kinetic energy. 
The first has no explicit spin interaction with the magnetic field while the second - the Pauli-Fierz 
model - has a ^/aa ■ B{x) term. The latter is much more delicate than the former (e.g., it requires 
bounds on a and Z, as well as the presence of a field energy, for stability, while the former needs no 
such restrictions), and our results in the second case are not as good as in the first. Still, they are 
meaningful. The third theory uses the Dirac operator and is relativistic (except for an ultraviolet 
cutoff in the A field). 

Clearly, it is essential to have a decent grasp of the binding energy, which is the truly physical 
quantity, in order to be able to proceed with a non-perturbative renormalization program. It is 
useful to recall that the binding energy for one electron and one nucleus of charge Z should be 
^mphysC^ a"^ Z"^ , where rriphys is the physical, renormalized electron mass. In the case of hydrogen, 
Z = 1, this is a comparatively tiny energy. (There is no charge renormalization in a theory, without 
pair production, but there is one in a theory with pair production, in which case a here must be 
replaced by the renormalized a.) 

Our purpose here is to find an upper bound to the binding energy in the two non-relativistic 
QED theories that have been extensively studied in the literature and to the relativistic theory in 
H. We believe it is the first time that such a rigorous estimate has been made. Note that the 
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binding energy (|T^) is always positive, by definition, and thus an upper bound corresponds to a 
lower bound on the energy of the fully interacting system relative to the energy of free electrons. 

If we equate the binding energy (which depends on the unrenormalized, or "bare" mass m) to 
the physical binding energy (which depends on rriphys) we obtain a (upper) bound to the amount of 



renormalization that is needed (see ( p^ ) and (27) for the first model). We believe these are the first 
rigorous renormalization estimates of their kind in a (admittedly simple-minded) quantized theory 
of electrons and photons. 

Let us note an interesting feature of our results about the mass. There are several ways to 
define the renormalized mass. The usual one is to look at the ground state energy of a free electron 
with fixed total momentum (electron + field) and to define l/2mphys to be the coeefficient of in 
the energy at p = 0. Another way is to set the binding energy equal to the physical binding energy, 
as is done here. 

• The latter definition has the property that for every value of A and of rriphyg there is a value of 
m that gives equality. (Our bound is unique, but the true answer is, conceivably, not unique.) 

• The former, usual definition very likely yields a solution for m only if A and l/m-phys are sufficiently 
small. (We always assume that m > in order that the notion of a ground state makes sense for 
the unrenormalized theory.) 

We cannot prove this last statement but it agrees with the prediction of perturbation theory and 
with classical electromagnetism. It also agrees with Van Kampen's exact solution of Kramer's 
"dipole approximation" model Q, which was the model that gave impetus to the renormalization 
program in QED. Instead of {p+A{x))'^ as in our first model, one takes {p+A{0))'^. In such a model 
the energy of free electrons goes as (7A^/^\/iV instead of CA^/^A, as in our model (see (pi])), 
but this is not the most significant point. For A = 1 one can compute the energy as a function of 
momentum p and finds that, with the first definition, mphyg = Cm + aA. Consequently, there is 
no solution for m if A or l/mphyg are large. 

The first Hamiltonian Hj^^ we consider is given by 

^^Hj^ = T + aVc + Hf, (1) 

where T is proportional to the kinetic energy of A^ electrons 

N 

r = Y,T, (2) 

with 

^^=^c^p^+^^^'^^^y ' (3) 

with pj = iVxj, and where m is the (unrenormalized) mass of the electron. The quantized, 
ultraviolet cutoff electromagnetic vector potential is 

A{^) = ^jzl («A(fc)e^'-^ + al{k)e-''-^) dk , (4) 



A=l 
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where A is the ultraviolet cutoff on the wave-numbers |A:|. The operators ax,a*^ satisfy the usual 
commutation relations 

[axik),al{q)] = d{k - q)6x,u , [axik), a^{q)] = 0, etc (5) 

and the vectors exik) are the two possible orthonormal polarization vectors perpendicular to k and 
to each other. The field energy is 

Hf=Y, I u,{k)al{k)ax{k)dk ; (6) 
A=l,2 -^K' 



the physical choice of uj is uj{k) = \k\, but our Theorems 2T and 3J are not restricted to this 
choice. No infrared cutoff is needed. 

Finally, there is the Coulomb potential. There are K nuclei with positive charges eZi, ...,eZK 
and with fixed locations Ri, Rk in M'^. In this model the nuclei will preferentially locate them- 
selves at those Rj that minimize the total energy, but these special locations are irrelevant for our 
theorem. 

■^■^ \xi — Rk\ ^ \xi — Xi\ ^ \Rk — Ri\ 

i=l k=l ' * l<i<j<N ' ^' l<k<l<K ' " " 

The Hilbert space is Hn = ^f=iL'^ {B? C'^) ^ T , where A denotes the antisymmetric tensor 
product (Fermi statistics), J- is the photon Fock space, q is the number of spin states for each 
electron (g = 2 in nature). The Hamiltonian (||) is bounded below by C{N + K) (even without 
the aid of Hj); this was first stated in |jl], p. 857] and p. 2], who noted that the stability proof in 



1 12] with ^ = extends to the A 7^ case by virtue of a diamagnetic inequality. In [U| a remark 



of J. Frohlich is presented, that this applies equally well to the quantized A field (0) because 



[A{x),A{y)] = for all x,y. We use and discuss this fact again in (17) below. 
The Pauli-Fierz Hamiltonian, , which is treated in Section ^, is 



^2 ^ 



hr 

hW = Hn + \/a Vo-,- • S(x,) with B (x) = curl A(x) . (8) 



2m 



The Pauli matrices Uj are the spin-^ operators for particle j. 

The third model is the relativistic no-pair model treated in Q, whose Hamiltonian is 



he 

Here, D(A) is the Dirac operator 



P+i^Y.D,{A)+aVc + Hf^P+ . 



l-Hj^' = P+{y]DM)+ «K + Hf}P^ . (9) 



771 C 

D{A):=a-{-iV + V^A{x)) + —(5 , (10) 

and P"*" is the projector onto the positive spectral subspace of D{A) for all of the electrons. 
Since the N Dirac operators commute with each other, this definition of P^ as a projector makes 
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sense. In other words, we start with the usual Hilbert space TCn and then replace it by the smaller 
physical Hilbert space 7i^^^ = P'^Hn- In W^^'' it is impossible to separate the spaces from 
the Fock space. 

While energy, being one component of a four-vector, is not a relativistically invariant quantity, 
it is true, nevertheless, that positive and negative energies of D{A) are relativistic concepts since 
they are invariant under Lorentz transformations that do not change the direction of time. We 
thank J-M. Graf for this remark and we thank J. Yngvason for noting that for this to be true it is 
essential that the joint spectrum of energy and momentum of D(A) lies in the light cone. We have 
not proved this, but it is plausibly true. 

2 Binding Energy Bound (simple version) 

In this section we analyze the binding energy for H^. 
We define Z = max{Zi, Zk} > and set 



K = -Z + 2.21599^/3^2/3 + 1.0307g^/3. 
We also define the (positive) free-electron ground state energy 



(11) 



Eo{N) = inf specFiv(O) > 



(12) 



(where H]\f{0) is the Hamiltonian (0) without the Vc term), the total ground state energy 



E{N) = inf specHN < Eo{N) 



(13) 



and the (positive) binding energy 



AE{N) = Eo{N) - E{N) . 



(14) 



THEOREM 2.1 (Binding energy for Hj\f). Assume one of two cases: 




Case A 



Case B. 



(15) 



[KaV2mc'^y^Eo{N)/N Case B. 

Proof. We use the known result for the stability of "relativistic" matter in the form given in [^, 
eqs. 2.9 and 5.2] (which improves some of the results in [|l3| ). 




N 




(17) 
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for any vector field A{x). (Note that although ( [T7| ) was proved for ordinary, non-quantized A fields, 
we are allowed to use it for our operator- valued field (Q) since the commutator ^(y)] = for 

all X, y, and hence there is a representation in which A{x) is an ordinary vector field. Of course, 
A{x) does not commute with Hj but that is immaterial.) 
By Schwarz's inequality 

Vc > -ny'2mc/hVNVT . (18) 

Since Hj > and since the square root is operator monotone, 

>T + Hf- Ka^j2mclh\fN^ >T + Hf- Kay'2mc/hVNy^T + Hf . (19) 

The final step is to note that the function f{x) = x — Kay^2mc/h\/lV ^Jx has its minimum at 
X = K^a^mcN /2h. Therefore, when x < n^a^mcN /2h we can say that f{x) is not less than the 
minimum of /(x), which is -\{KafmcN /h. This is Case A of (|l5|). Therefore < Eo{N) + 
^{Ka)'^mc^N < {Ka)'^mc^N, which is Case A of (0). On the other hand, if x > n^a^mcN /2h we 
can say that f{x) is monotone increasing in x. Since x = T + Hj > Eq{N)/11c > n^a^mcN /2h, 
the infimum of the spectrum of the right side of ( p!9| ) is given by f{EQ{N)/hc), which is Case B of 
(0). □ 

To apply this Theorem 2.1 we must have a decent estimate of Eq{N). Let us consider the 
physical case u;{k) = \k\ and let us define (with Ac = ^c/mphys = physical Compton wavelength) 

7^=^ = ^AAc, (20) 

which is the ratio of the cutoff photon energy to the self energy that an electron would have in a 
relativistic theory. A bound on Eq{N) in this case is provided in where a proof is announced 
and outlined that for fermions there are constants Ci,C2 (depending on q) such that (for large A 
and fixed a) 

1/2^3/2^ < Eo{N) < Ci m(?o?l'^Tl^'^I^N (21) 

The exact exponent is still not known but we lean towards 12/7. In any case, it differs from the 
perturbation theory value 2. Fermions are most important here because one can show Q that 

1/2^3/2^1/2 < < 2/7^12/7^5/7 foj, 

bosons, and this would be useless for 
our purposes. Unfortunately, the bounds in (21) do not imply that Eq(N) is strictly linear in N, as 



one would hope. We also note that ( pij ) holds even if the Coulomb repulsion among the electrons 
is omitted. 

If TZ is small then we are in Case A. This will surely be so if K^a 

12/7 > 2C72 7^l2/7. Let us note 

that K, < 5.67Z for Z > 1 and q = 2. Then, 



< (32.2)a^Z^mc^ Case A , (22) 
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which compares not unfavorably with the hydrogenic value ehydrogonic(-^) = ■2^^ct^^physC^/2, where 
mphys is the physical electron mass. As A increases the bare mass m should decrease. If we set 
AE{N)/N = ehydrogenic(-2^)) inequality (22) tells us that the required bare mass cannot be too small, 
namely 

m > mphys/64.4 Case A . (23) 

We turn now to Case B, which surely holds if 

K2a3/2 < 2C7l7^3/2 . (24) 

With the help of (0) we can conclude that 

< Ka^^n^^mc^ Case B . (25) 

To understand Case B further, let us use n < 5.67Z and note that (^) becomes 

AE(N) 5.67 C^^ R/7 m ^ 
^^7^<^47y2i3/V/' if 5.67Z</., 26 

A'ehydrogenic(^) Z (^4/7 fUpi^ys 

with fi = a-3/4^/2c77e3/4^ 

To satisfy the condition in ( p6| ) for all Z < 92 we can take fi = 5.67 • 92 = 522 or 7^ = 
(522)''/3(2Ci)-2/3q, = 30.7(2Ci)-2/3 (with a = 1/137)). This means that we fix A in units of the 
bare Compton wavelength hc/m. Assuming that we choose Ci to be not too large (which can always 
be done since Ci refers to a lower bound in (21)), this allows for a sizeable value of the cutoff A 
(see (pO|) ). Now let us set the left side of (|2^ ) equal to 1, in order to make contact with experiment. 
We then find (since Z > 1) that 

m > mphys Cf^'^C^^^'^/lS, 800 Both Cases . (27) 

Alternatively, we may measure A in terms of the physical Compton wavelength Ac. That the 
bare mass cannot be too small can be seen as follows. Consider the following inequality, which is 
related to (pi) 



(5.67 • 92)^0^/2 < 2(7l7^3/2 . (28) 
If this inequality fails then we have the bound 

> mphysAAc(2C7i)2/3(5.67 • 92)^^/3^-^ = 1.74 • 10-'^(2Ci)2/3AAcmphys (29) 



m 



with a = 1/137. 

On the other hand if ( [2^ ) holds then ( p4[ ) holds for all Z < 92 (since k < 5.67 Z). Then we are 
in Case B and if we express the right of ( p6|) entirely in terms of AAc we find (for all Z < 92) 



m > (2-5.67y2^)-V(92)^(AAc)~^mphys (30) 
= 3.0 X 10-^^2^'/' (AAc) -'mphys • (31) 
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3 Binding Energy Bound (Pauli-Fierz version) 

In this section we analyze the binding energy for in (^). The ground state energy and binding 
energy are defined as before in (0), (|l|), (|T|), but with H^^ , and we do not encumber the notation 
with a superscript PF. As far as constants are concerned, the following theorem is not the best 
possible one, but it is presented this way for simplicity. In particular, we do not have to assume 



that a is bounded — as we do in the hypothesis of Theorem 3.1. 

Some constants have to be defined. The maximum nuclear charge Z is defined as before and 
we then define O = max{Z, 20.6}. (Note, for later use, that 20.6 = 64.5/7r.) We also define 
C = (0.060) (Stt) (3/4) -3/2 = 2.322 and 



aK 



1/4 



r = 9.65 ^—j A . (32) 
We also define the operator hcT' to be the total Pauli-Fierz kinetic energy, namely, 

T' = j:Tr = T + ^y^E • Bi.,) > , (33) 
3=1 i=i 

where 

^^''^ = ^ci^P' + ^^(x,))' + V^<rr B{x,)} = ^ {a, . {p, + V^A{x,))Y . (34) 

THEOREM 3.1 (Binding energy for H^^). Assume that lirQ^a^ < 1 and assume one of two 
cases: 

' < 2(1 - 27rJl^a2)-2(7rJla)2 mc^N, Case A 

Eo{N) { (35) 
> 2(1 - 27rf]Ca2)-2(^j^^)2 ^^'^jy^ Case B. 



Then the binding energy per electron satisfies 

2(7rJla)2(l - 27rJlCa2)-2(2 _ 27r0^a2) mc^ + 2TrnaThc Case A 
27rJ]^a2 Eo{N)/N + 2-KnaV2mc'^y^Eo{N)/N + 2TrnaThc Case B. 



AE{N) ^ 



N 

Proof. The strategy is the same as in Theorem 2T. An analogue of (|T7| ) is provided by @. 



(36) 



1 ^ 

> - E ■ {Pj + V^A{x,)) I - ^aHf - TN , (37) 

The derivation of (^^ from |^] will be explained at the end of the proof. For the moment let 
us continue with the proof of the theorem. The analogue of (O) is then 



> -\/2mclh\fN\fT - ^aHf - TN . (38) 
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Consequently, since T' > 0, 



(39) 



The rest follows as in the proof of Theorem 2.1 



It remains to show how (37) arises from ref. Q. We consider a Hamiltonian, Hi\f, similar to 



H'^ in but with some auxiliary parameters. 



N 



Hn = Y1 ki • {Pj + VaA{xj)) I + a'Vc + -fHf 



(40) 



with a' , 7 > 0. Note that Hjsf has the dimension of length"^ and not an energy. Note also that 
the a appearing in pj + ^/aA{xj) is the true a and not a' . 

We shall prove the analogue of Theorem 2.2 of |^] with Hn in place of H'^ and with an 
appropriate substitute for the lower bound stated in Theorem 2.2. Theorem 3.1 of is unchanged, 
and we turn to Sect. 4. We take e = (it can be taken to be zero even if m 7^ 0). All the equations 
in this section remain true if we replace ku by Ka' and k^o^ by K^(a')^. Thus, in the lower bound 
for the operator H2 in Sect. 4 of Q, there is a in the numerator and a' in the denominator. 

The next step is to use the inequality in Example 1 of Appendix B of Q to bound J appearing 
in the lower bound for H2, but we have to remember that we have and not Hf. Thus, the 
analogous conditions on the parameters are k < irQ, Ka' < 1 and 7 > 87r(0.060)a (l - K^ia'f) ' . 

We make the choices 

At = 7rO, a' = (27rl7)-\ 7 = ^0. (41) 

The lower bound to Hjq is as in Q. We find that the corresponding number C2 is bounded by 

6+ (a72)(V2Z + 2.3)2 



'k 



(27/2vr)7 



< 



N 
K 



137r 
277 



(42) 



since (\/2Z + 2.3) V27rQ < 1. 

The analogous lower bound for the operator H^, is then 



A^^ > "A 



"18A7' 


cl>-{ 


vr 





aK 



1/4 



with V = {i/nf^ (2/3)(39)3/4 = 9.65. 

This completes the derivation of (|37|) from 



(43) 



□ 



Using Theorem 3.1 we could proceed to derive more explicit bounds for the binding energy and 



mphys - as in the discussion after 2J, but we leave this task to the interested reader. The only 
needed information is the analogue of (p]) for the Pauli-Fierz operator. As announced in 



a 



1 + a2/3 



(44) 
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4 Binding Energy Bound (relativistic version) 

We consider the Hamiltonian in (P). Our results here are very crude and we state them 
mainly to point out that realistic results on the binding energy could be obtained if one were 
able to improve the estimates of various constants. In the present situation we do not have any 
bounds on Eq{N) (other than the simple one Eq{N) > mc^N), but we expect something like 
Eq[N)/N ~ mc^ + (const. )/icA for not too small A. 

We set T" = P+ J2Zi I A(vl)|P+. With P+D{A)P+ in place of \cjj ■ {pj + ,/^A{xj)) \, @) is 
valid for P~^VcP~^ — in the same way that the inequalities of Theorems 2.1 and 2.2 of are valid 
with the same constants. We find that 



—H',^^ > -2-KVLaT" - - 2-nVLaTN + T" + Hf 

he 

= (1 - InQa) T" + (l - 211^^0'^) Hf - iTrQaTN (45) 
> (1 - InnCa"^) [T" + Hf) - 2-K9.aTN 



since < 1. Therefore, 



^<2.na^ + 2.naTfu:. (46) 



We note that — apart from the unnaturally large constant — the binding energy appears to be 
bounded by a times the self-energy. We also note that (45) and ( |4^ ) can be improved a little by 



using the free parameter < e < 1 that appears in Sect. 4 of Q; we have taken e = here, as we 
did in Sect. ^. 
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